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Abstract 

Based on a field theoretically inspired model of light-cone wave functions, we derive valence-like 

generalized parton distributions and their double distributions from the wave function overlap 

in the parton number conserved s-channel. The parton number changing contributions in the t- 

channel are restored from duality. In our construction constraints of positivity and polynomiality 

O ■ are simultaneously satisfied and it also implies a model dependent relation between generalized 

] parton distributions and transverse momentum dependent parton distribution functions. The 
O ■ 

• • ■ model predicts that the t-behavior of resulting hadronic amplitudes depends on the Bjorken vari- 



able XBj . We also propose an improved ansatz for double distributions that embeds this property. 
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i. Since the discovery that the proton is a composed system, enormous amount of efforts have 
been made in order to reveal its contents and understand the dynamics of its constituents. Var- 
ious frameworks and models have been simultaneously proposed; often based on the quantum 
mechanical concept of the proton wave function. Factorization theorems, derived in the frame- 
work of perturbative Quantum Chromodynamics (QCD), are the basis to relate non-perturbative 
parton distribution functions (PDFs), distribution amplitudes, and generalized parton distribu- 
tions (GPDs) with experimental observables. It is obvious that these non-perturbative quantities 
are somehow related to the proton wave function. However, quantifying this relation starts with 
the conceptional difficulty of defining the bound state problem in the relativistic quantum field the- 
ory, and it ends up with the difficulty in matching the evaluated (generalized) parton distributions 
with those that are defined implicitly in the perturbative factorization approach. 

The idea to write down the proton wave function in terms of the partonic degrees of freedom was 
spelled out in the early days [1, 2, 3]. Thinking of the proton as a bunch of partons that move nearly 
on the light-cone, e.g., specified by = (1,0,0, —1), allows more easily to establish the desired 
link, see Refs. [4, 5] and references therein. This leads to the concept of light-cone (LC) wave 
functions ^/^^'^(Xj, k_|_j, Aj). They are the probability amplitudes for their corresponding n-parton 
states \n,pf ,p±i, Xi), which build up the proton state with LC hehcity S = {+1/2(1), —1/2(|)}: 

r " 1 

\P,S = {T, i}) = J2 ^'"^l" "f^n^Xi, k^. A,) J] -y= |n, X,P+, XiP^ + k^„ A,) , (1) 

where we used a shorthand for the n-parton phase space: 

li,A,l = fl'-^^16.H(^l-±x}jS^^^(±i.,}j . (2) 

The LC wave functions depend on the longitudinal momentum fractions Xj = k^/P~^ (the plus 
component of a four- vector is V'^ = + = n ■ V) , the transverse momenta k^j, and the 
LC helicities Aj. They are determined from the eigenvalue problem for the LC Hamiltonian P~: 

p-\p^S) = —\P,S) , with P- = P°-P3^ P+ = P° + P^ P± = 0, (3) 

which can be derived from the QCD Lagrangian. However, in practice the QCD dynamics is not 
well understood and it remains very challenging to develop this concept to a stage at which it can be 
used for quantitative evaluations of physical observables or parton distributions [5]. It is common 
to pin down LC wave functions by comparing their resulting model predictions with experimental 
observations. Certainly, so far this concept is extensively elaborated, it connects different non- 
perturbative quantities, like PDFs, transverse momentum dependent parton distribution functions 
(TMDs), distribution amplitudes, form factors, and GPDs to a more universal object. 
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Our study is mainly devoted to GPDs, which are accessible from hard-exclusive reactions, 
e.g., electroproduction of mesons and photon. They arise from the non-diagonal overlap of LC 
wave functions, and therefore contain a maximum of information about the proton wave function, 
compared to other non-perturbative quantities. For instance, they are reducible to elastic form 
factors and PDFs. Field theoretically they are defined as off-diagonal matrix elements of two field 
operators that live on the light cone [6, 7, 8]. Describing the initial and final proton states, with 
given momenta Pi and P2 = Pi — A and LC helicities 5*1 and S2, respectively, in terms of the 
LC wave functions (1), one can straightforwardly derive wave function overlap representations of 
GPDs for the partonic s-channel exchange [9, 10, 11]. In this partonic process the number of 
partons is conserved and the momentum fraction x of the struck quark is larger than the skewness 
parameter rj = {Pi — P2) / {Pi + P2) > (up to a minus sign we use the variable conventions of 
[6]). In this outer region of x, the GPDs H and E (Ji's conventions [12]) read 

2 - C \- \- / T2-n 



H-j^E]{x>r],r],t)= (4) 



n A 

Ai-tA2 
2M 



/W'"" / [dX d'k^l 5{X - XO^I*(X:, k'^,, XMiX^, k^„ A,) 



E{x>7],V,t)= (5) 
V^YY J [dX rf'kx]„ S{X - Xi)i;l*{X:, k'^,, \M{X,, k^„ A,), 

n Xi 

where ( = 2ri/{l + 7]), X = {x + ri)/{l + t]), and the momenta of the outgoing partons are 

^ ]^ _ 

X[ = — — , k'j_j^ = k_|_i ^^-L struck quark, (6) 

X X 
X[ = - — ^— , k'_i_j = k_Li + - — ~^^-L ^'^^ spectators i = 2, . . . ,n. (7) 

Anti-quark GPDs are analogously defined with a negative momentum fraction x < —r]. The 
central region, i.e., —rj < x < rj, arises from the t-channel process in which the parton number 
changes from n + 2 to n [10, 11]. Viewing this as a mesonic-like t-channel exchange makes contact 
to Regge phenomenology [13, 14, 15]. Note that positivity constraints, in its most general form 
[16], should be satisfied in the overlap representations [17], if they are not spoiled by subtraction 
procedures. Indeed, this can be easily shown for a two-body LC wave function, as used below. 

Let us also remind of the constraints of Lorentz covariance for (quark) GPD form factors. They 
are not invariant under general Lorentz transformations, however, they are built by a series of local 
twist-two operator matrix elements, belonging to irreducible representations, that are labelled by 
the spin J > 1. It turns out that Mellin moments of GPDs H + E and E with the weight x'^^^ 
are polynomials in 77 of the order J — 1 and J, respectively. Time reversal invariance combined 
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with hermiticity requires that these polynomials are even [12]. These properties are manifestly 
implemented in the double distribution (DD) representation of GPDs [6, 7]: 

where the DDs e and h are symmetric in z. The DD representation for E is not uniquely defined 
[18, 19]; Eq. (8) shows the one which naturally occurs in our model studies, see below Eq. (20). 

As explained above, the partonic interpretation of GPDs separates the support in the central 
(|x| < 77) and outer (77 < |x|) regions. Residual Lorentz covariance, ensuring the polynomiality of 
moments, requires that both regions are tied to each other and that the functional form of GPD 
is constrained, e.g., in the outer region the GPD E is given by the integral representation: 

„£±22 , 

/ cly 

E{x > r],r],t) = {1 - x) — e{y, [x - y) /r],t) . (9) 

l-r] 

We consider both regions as dual to each other, i.e., knowing a GPD in one region allows to 
restore it in the other. Hence, a GPD can be entirely evaluated from the parton number con- 
served s-channel overlap of LC wave functions. A constructive, however, unwieldy method for the 
restoration of the central region is based on Mellin moments and its inverse transformation [14]. 

In the following we utilize two-particle LC wave functions. They serve us to describe the proton 
contents by a constituent quark and scalar diquark, where the latter plays the role of an collective 
spectator. Numerous investigations in this spirit, even much more advanced ones with specific 
emphases, e.g., Refs. [20, 21, 9, 22, 23, 24, 25, 26], have been made in the past. The new aspect in 
our study is that we take care of Lorentz constraints for the LC wave functions, cf. Refs. [27, 28]. 
This allows us to evaluate^ DDs from the parton number conserving overlap representations (4) 
and (5). We outline the generic features of such model, illuminate their t-dependence, point out 
its restrictions, and overcome them by hand, yielding improved DD and GPD ansatze. 

ii. The functional form of LC wave functions is dictated by the underlying Lorentz symmetry, 
i.e., the longitudinal and transversal variables are tied to each other in a certain but not apparent 
manner. Hence, writing down LC wave functions by hand usually results in a violation of the GPD 
polynomiality property. Note that this failure can not be fixed by taking into account the particle 
number changing processes. The guidance for an appropriate model comes from a perturbative 
calculation in lowest order [30]. We employ the Yukawa theory and have the LC wave functions 
for four helicities, 

^/jly.iX, k^) = ^U^,{X, kx) = (m + ^) cpix, kx) , (10) 

V^|,/,(X, k^) = -^^{X, kx) , £,/,(X, k^) = ^^(^, k^) , (11) 

^Of course, this task is straightforwardly done in a Lorentz covariant formalism [29]. 
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in terms of a scalar function ip(X,'k±). This scalar function arises from the spectator propagator 
in a triangle Feynman diagram [30, 31] and so the underlying Lorentz symmetry is respected. We 
generalize (p{X, k^) by an adjustment of its power behavior p: 



gM^P 



(12) 



X 1-X 

where M, A and m are the proton, spectator, and quark masses, respectively. Note that the factor 
X'P takes care of the proper Lorentz behavior and that the Yukawa theory result is for p = 0. 
The GPDs are now evaluated in the outer region from the overlap representations (4) and (5), 

^e) ix>r],7],t) = (13) 



H - 



7] 



2-C 



2M 



2^1^ J 167r3 L 
E{x > T], r], t) = 



^Pl\/,ix', k'^)ijl,/,{x, k^) + i^%{x', K)i^l,/,{x, kj 



(14) 



— ^ijl\/,{x\ ki)v;ii/2(x, k^) + ^%,ix\ kl)V;ii/2(X, k^) 



for the two body wave functions (10)-(12). For p > the k_L integrals are finite, while for p = 
the representation (13) suffers from an ultraviolet divergence. We find for the GPD E: 



E{x > T], rj, t) 



2r(2p+l) 



da 



[(1 - X)(l - X')r' [aar (f + X - a(l - X%) 

Afmin-fl2P+ 



(47r)2 T{p + 1)2 7o [f{X\a) + f{X'\a) + aa{l - X){1 - X')'-^] 



where U 



—C'^M'^I (1 — C), a = 1 — a, and the mass terms are collected in 



f{X\a) 



a 



;i-x) 



m 



\2 

X—--X(l 
M2 ^ 



X) 



T, (15) 



(16) 



The result for H has a similar structure and will not be displayed for shortness. 

Since our model respects the underlying Lorentz symmetry, there must be now a possibility to 
transform the overlap result (15) into the form of the DD representation (9). In fact, this can be 
simply achieved by a linear variable transformation of the integration parameter 

1 1 — r] ( ^ X — y 



a 



21 -X 



y 



V 



(17) 



and removing the residual skewness dependence by using x = y + -qz. We arrive at 



E{x > r], r], t) = (1 

where the DD is given by 

eiy,z,t) = N 



i+'7 dy 

l-ri 



x-y_ 
V 



[(l-y) 



j^ + y^-yi^- y) - ((i - y? - 



(19) 



Here we absorbed several factors, including (7^, in the normalization constant A^. From the DD 
(19) and Eq. (8), we find E for the central region, arising from parton number changing processes: 



i+'i ay I X — y 



E{-r]<x<r],r],t) = {l-x) / —eh/, "-,t] . (20) 

Jo V \ V J 

The evaluation of the GPD H in terms of the DD h goes along the same line. The combination 
{H + E){x, 7], t) can be written in the form of the DD representation (8) with 

.l-2p {{l-yf-z^)P 



{h + e){y,z,t) = N- 



+N- 



[(1 - y)^ + y^- -y)- ((i - yf - ^')mA 

^-^ + y + y^ + {2-y)f^] {{l-yf-zy 



2p 



2[il-y)^ + y^- yil - y) - ((1 - yf - z^)^] 



2p+l 



(21) 



Finally, the overall constant is fixed by the normalization condition: 

dx H{x,ri,t = 0) = dy dz {h + y e){y, z,t = 0) = 1 . (22) 
-1 Jo J-l+y 

iii. We employ now our simple-minded model to evaluate form factors, TMDs, PDFs, and GPDs 
as net contributions of valence-like quarks in the proton, i.e., we deal with the isospin combination 

H{x,ri,t) = Hu,^X^,r],t) - Hd^^^{x,ri,t) , E{x,r],t) = E„^^,(x, r/, t) - Ed^^^{x,r],t) . (23) 

We have only three model parameters, i.e., power behavior (p) of the LC wave function, spectator 
(A) and quark (m) masses. Our goal in doing so is to explore this model and find its restrictions. 

We start with the electromagnetic form factors, which are obtained from the Drell-Yan-West 
formula or, equivalently, by integrating out the momentum fraction dependence in GPDs: 



Fi{t) = j dxH{x,ri,t), F2{t) = j dx E{ 



x,r/,t). (24) 

The Dirac form factor is normalized by -Fi(t = 0) = 1, cf. Eq. (22). The asymptotic drop off for 
large —t is estimated to be l/(— t)^, according to field theoretically inspired [2] and phenomenolog- 
ical [3] model counting rules and the perturbative analysis [32]. This suggest to set p = 1. We use 
the charge radius squared = QdFi{t) / dt\t=Q and the anomalous magnetic moment k = F2{t = 0) 
of the proton to pin down the remaining two parameters. The following plausible mass values 
yield an agreement with experimental measurements, cf. Ref. [20]: 

p=l, A = 0.75GeV, m = 0.45GeV i? = 0.76fm, k = 1.78. (25) 

The form factor Fi{t) fairly agrees with experimental data as shown in Fig. 1(a) and behaves 
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Figure 1: The form factor Fi{t) (a) and the ratio {—t)F2{t) / KFi{t) with a cut-off for the endpoint 
singularities (b). Data are taken from Refs. [33, 35, 36]. 

for —t ^ 100 GeV^ as t)^, however, in the large —t asymptotics it drops faster. We would 
conclude, in agreement with aforementioned counting rules, that F2{t) decreases as l/(— t)^. This 
follows from counting the powers of in Eq. (14), arising from k^-dependence of wave functions 
(10-12), and the vanishing of the overlap integral for E if spherical symmetry is restored in the 
asymptotic A^^ oo. However, we observe that the limit A^ — s> oo can not be taken before the 
k^ integration, since divergences appear and render the integral to be infinite. This behavior shows 
up also in the DD (19), where the limit —t^oo causes end-point singularities ai z = ±(1 — y). 
Hence, we effectively find a F2{t) ~ l/(— t)^ behavior and thus naive power counting fails. We 
remind that the assumptions for the counting rules were carefully spelled out [2, 3, 32]. 

End-point singularities arise also in the perturbative evaluation of F2 and their regularization 
yields a logarithmical modification of the l/(— t)^ scaling [37, 38]. The experimental measure- 
ments indicate a F2{t)/Fi(t) oc l/-\/^ scaling [34, 35, 36], which might be also interpreted as 
a logarithmical modified F2(t)/Fi(t) oc l/(— t) scaling. Interestingly, the anomalous break down 
of the naive power counting in our scalar diquark model leads to a F2(t)/Fi{t) oc const, scal- 
ing. Applying the recipe of end-point regularization to our result, i.e., imposing the constraint 
0.33 < 1 — y — \z\, we easily can fit the experimental data as demonstrated in Fig. 1(b). This 
exercise should not be considered as a serious attempt to explain experimental data, rather as a 
demonstration showing that fitting can be done easily. One might wonder whatever this failure 
of the scalar diquark model is related to the oversimplified spin coupling and could be cured by 
inclusion of an axial-vector diquark, or it might simply reflect a wrong implementation for the 
large k^ behavior of the LC wave functions. 

Our model also predicts TMDs and their k^-integrated PDFs outcome. Note that the TMD 
concept has numerous fundamental issues, see, e.g., Ref. [39]. We might define here a unpolarized 
valence-like TMD in terms of the LC wave functions (10-12) overlap, where C = ^ind t = 0: 



q{x, kj 



g^M^P [{xM + 



m] 



N2p+1 



[k^ + (1 - x)m2 + - x{l - x)M2] 



2p+2 ■ 



(26) 
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At large they fall off with l/(k5^)^^+^ and are suppressed in the large x region by (1 — x)^^"^^. 
Our model is similar in spirit to the spectator model utilized in Ref. [40] for the evaluation of 
kx-(un)integrated parton densities and fragmentation functions. There the fermionic propagators 
are replaced by quark-diquark form factors with a cut-off mass A, while in our case they are take 
to be on-shell. Replacing m by A, we find that Eq. (26) is identical with the result (80) in Ref. [40]. 
The valence-like PDF is obtained from (26) and is related to the GPD H: 

gix) ^ Hix,, = 0,t = 0) = Tx2^n" ITJp^^ ' ^''^ 

2p{2p + 1) (xX^ + {1- x) {m? - M^x))^^ 

For p = 1 we obtain the generic behavior of parton densities at large x. Note that the form factor 

Fi{t) falls at large —t with (— 1)~^~^. Hence, setting p = n — 1 we confirm the known counting 

rules. The momentum fraction carried by the valence quark combination u — d, 

(x) = [ dx X q{x) , q = qu^^^ - qa^^^ , (28) 
Jo 

is (x) ~ 0.26 with our parameter specification (25). It nearly agrees with the value in the Gliick, 
Reya, and Vogt parameterization [41], given at a low input scale /Iq = 0.4 GeV^ (to perturbative 
next-to-leading order accuracy): {x)u^^i — {x)d^^i ~ 0.24. This amazing coincidence of momentum 
fractions should be considered rather as an accident. For x going to zero, the PDF (27) approaches 
a constant. This behavior we consider as an unrealistic feature of the diquark model. The small 
x-region, i.e., the high-energy limit might be understood in the Regge picture as an exchange 
of mesons in the t— channel which leads to the expectation that the valence-like parton densities 
behave as x~"'^^\ where a(0) is the (effective) intercept of the meson trajectory. From the s-channel 
view, which we take, the true small x-behavior arises by summing up all Fock state components. 

We come now to the GPDs. First, we comment on our duality assumption between the partonic 
s-channel and mesonic like t-channel exchange. We could have added to the central region a term 

D{x/ri,t) =e{\x\ <\ri\)d(-,t] , with d{l,t) = d{-l,t) = (29) 



that entirely lives in the central region, vanishes at the cross over point and is anti-symmetric in 
X. An explicit GPD evaluation from the parton number changing processes confirms that such 
a term is absent and so the underlying field theoretical model respects duality. However, such a 
D-term is needed to cure the common DD representation for the GPD E or H [42]: 



^|(x,?7,t) = y dy J dz 6{x -y - r]z) <^^j {y, z,t) ± D{x/r],t). (30) 

Our finings suggest that both e^''^'^{y, z,t) and d{x/ri,t) have a cross talk. Indeed, they are two 
different projections of our DD (19); the D-term is extracted from t] ^ oo with fixed x/f] [43]: 

.(., t) ^ X r\y iV(f + .)((l-.)--xY 

[(1 -y)jr^ + y^-y{^-y)- ((i - yf - 
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Figure 2: The GPD H{x,r],t) in the spectator model (8,19,21) (a) and factorized VGG ansatz (30,33) 
(b) with the parameter values in (25). The momentum transfer squared is set to t = imin — 0.25 GeV^. 

Interestingly, the coefficients in its Gegenbauer expansion [44] are rather small and sign alternating: 

d{x, t = 0)^(1- x^) [0.345 Cf/^(x) - 0.163 C'l^\x) + 0.026 C'l^^x) + ■■■]. (32) 

The first coefficient is in astonishing agreement with the estimate of the chiral soliton model [45] , 
given as 0.33 at a intrinsic scale of /ip ~ 0.36 GeV^. Note that the flavor singlet quark D-term in 
the chiral soliton model is predicted to be negative and large in magnitude [44, 46, 45]. 

In confronting GPDs with experimental data, one usually factorizes^ the t-dependence in the 
DD ansatz [47, 48, 44] (VGG refers to the popular code of Vanderhaeghen, Guichon, and Guidal): 

h-^^^iy^z^t) = F,(t)^nVGC(,/(i _ ^)) , nVG«(.) = -^^^Al - z^r. (33) 
l-y V7rl(o+lj 

The even proffie function n^^*^(z) is chosen to be convex and normalized to J^^dzlV^'^'^{z) = 1. 

Let us compare the spectator model (8,19,21) with the t-factorized VGG ansatz (30,33). The 
latter we make from the PDF (27) and D-term (32), multiplied with the form factor (24). For the 
proffie function we take n^*^^(z) = (3/4)(l — z^), which also appears in our DD (21) with p = 1 at 
t = 0, cf. Ref. [47], and employ in both models the parameters (25). In Fig. 2 we show the shape 
of GPDs versus x and 77, where the momentum transfer squared is set to t = tmin(0 ~ 0.25. The 
differences between the spectator model (a) and the t-factorized (b) GPDs are clearly visible at 
larger values of 77, where the later ones are broader in their x distribution and have a smaller value 
[see also below Fig. 3(c)] on the trajectory x = r], compared to the former ones. The important 
difference between the two models is in their t-dependence. In Fig. (3) we display in panel (a) 
the t-dependence of the ratio H{ri,ri,t)/H{'r],ri,0) on the trajectory x = 77 for various values of 
77. Obviously, the t-dependence is flattering out for larger values of t], while in the small t] region 
it even becomes steeper, compared to the t-factorized GPD ansatz (dashed line). Analogously, 



^Such an ansatz might be improved by replacing the effective Reggc intercept in the quark densities by the 
dependent Regge trajectory, which wiU only partly remove the factorized t-dependence. 
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Figure 3: The t-dependence of H{x,7],t)/ H{x,r],t = 0) at x = r/ in the t-factorized ansatz (30,33) (thin 
dotted) and the spectator model within t] = {O.l(dotted), 0.3(dashed), 0.5(dashdotted), 0.7(sohd)} (a) 
and for the Mehin moments n = {O(dotted), l(dashed), 2(dashdotted), 3(sohd)} at 77 = in the spectator 
model (b). Imaginary (c) and real (d) part of the amplitude (34) versus ^, arising from the spectator 
model (solid) and factorized ansatz (dashed), where the dotted curve shows the real part without D-term. 

we find for the ratio J^dxx"'H{x,0,t)/ J^dx x^H{x, 0,0) of Mellin moments at 77 = that the 
t-dependence becomes flatter with increasing spin n as shown in panel (b). Such a behavior was 
also seen in lattice calculations, see Ref. [49] and references therein. 

As in the case of PDFs, the small x-behavior of GPDs in a spectator model should be considered 
as unrealistic. Since at x = 77 the momentum fraction X' vanishes in the overlap representation 
(13,14), we realize that on this trajectory realistic GPDs can be obtained only if the small X 
behavior of the wave functions is understood. This also means that one has to sum up all Fock 
state components, see discussion in Refs. [9, 28]. Note that already the evolution to leading order 
in the flavor singlet sector knows about the small x behavior, however, in the non-singlet sector 
one has to perform a resummation of t-channel exchanges, perhaps, along the line as suggested in 
Ref. [50]. For larger values of x on x = rj one is forced to understand at the same time the large 
and small X behavior of the wave functions and their interference. 

Having this warning in mind, we have now a look at the resulting amplitude, which appears 
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in the hard exclusive photon or electroproduction to leading order of perturbative QCD: 



The imaginary and real part of the amplitude versus ^ are shown in Fig. 3 (c) and (d), respectively, 
where we set t = tmin — 0.25. Note that the imaginary part is given by 'kH{x = C,, C,, t) and provides 
us the GPD shape on the trajectory x = C,- The differences between the spectator (solid curve) 
and t-factorized (dashed curve) models are obvious. Both the imaginary and real part for the 
spectator model is (much) more enhanced at large C,- For ^ 1 the imaginary part (3(c), solid) 
behaves as (1 — ^), while for the t-factorized GPD ansatz an additional suppression factor appear: 
(1 — ^)(— 1)~^ ^ (1 ^ 0(~^min)~^ ~ (1 — On the other hand, we observe that at smaller 
values of ^ the imaginary part of the spectator model is getting smaller. We stress, however, that 
small ^-physics is not implemented in the model. The sign change of the real part in panel (d), 
somewhere in the valence quark region, is a feature that is also observed in other valence-like GPD 
ansatze. The position of the zero, however, is floating and strongly depends on the chosen ansatz. 
At smaller values of ^ the real parts of both models approach to each other. The dotted curve 
shows the t-factorized ansatz without D-term, leading only to a slight constant shift. 

Concerning the results of our spectator model, we conclude that the t-dependence of cross 
sections for deeply electroproduction of mesons in the large Xbj region is getting flatter and that the 
cross section could be (much) larger than the prediction from a t-factorized GPD ansatz, cf. Figs. 3 
(c) and (d). This is consistent with preliminary measurements of the e~p{Pi) — > e~p{P2)p^^^ cross 
section, where the t-slope decreases if larger values of xbj and are approached [51]. We find, for 
instance, that the slope for the amplitude square |7i(.^,t)p, parameterized as e^*, decreases from 
b ^ 3.5 GeV~2 at xbj = 0.2 and = 2 GeV^ to b ^ 1.5 GeV"^ at xb^ = 0.6 and = 5 GeVl 

iv. In this paper we derived from the overlap representation of LC wave functions the valence-like 
GPDs and their relatives: proton form factors, TMDs, and PDFs. We generalized the LC wave 
functions from a scalar diquark spectator model in such a way that the non-manifest Lorentz 
behavior of the LC wave functions is respected. This is the key to obtain the DDs from the 
overlap representation in the partonic s-channel and then to restore the full GPD support. 

Our model fairly describes the Dirac form factor Fi(t) and the proton anomalous magnetic 
moment comes out correctly by a natural choice for the constituent quark and diquark masses. 
However, the model fails in the description of the t-dependence for the Pauli form factor F2{t). 
Namely, its naive power counting for the large t-behavior is spoiled by end-point singularities that 
appears at the branch point —t = oo. Unpolarized valence-like TMDs and PDFs, also obtained 
in Ref. [40], have the expected generic behavior at large momentum fraction x and the average 
value {x) fairly agrees with the result of Ref. [41], given at a low input scale. 




(34) 
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The GPD models satisfy by construction the positivity and polynomiahty constraints. The 
found DD representations (8) naturally completes the polynomiahty condition and avoids so a 
'misleading' D-teim phenomenology. Another important characteristic property of the model is 
that the t-dependence of GPDs, resulting from the DDs (21) and (19), is washed out at large x. 
This behavior is rooted in the fact that the variables t and x are intrinsically correlated because the 
transverse and longitudinal degrees of freedom in the wave functions are tied by hidden Lorentz 
symmetry. This flattering of t-dependence also appears on the trajectory x = rj and therefore it 
can be confronted with experimental measurements. Note that the t-behavior of GPDs is only 
tested in lattice calculations for 77 = 0. From the overlap representation it is clear that the t- 
dependence of GPDs and the k^-dependence of TMDs are closely related to each other, since 
both arise from the k^-dependence of wave functions, a recent discussion is given in Ref. [52]. 

It is in the nature of a spectator model that it fails to describe the small momentum fraction 
behavior, which arises from the summation over all partonic Fock state components. Hence, there 
is a potential problem for the predicting power of such models for GPDs on the trajectory x = t], 
accessible in experiments. Here a GPD arises from the interference of the LC wave function at the 
extreme limit of vanishing momentum fraction with that of a non-vanishing momentum fraction, 
controlled by the skewness parameter. Therefore, even the limit x = 77 — > 1 is rather intricate [53] . 

It remains a challenging task to construct GPDs that satisfy all theoretical constraints and 
are flexible enough for a 'global' fit of experimental data. For the time being, we might suggest 
to adopt the t-dependence in the ansatz for the DD and improve its failure at small y, i.e, small 
X for the resulting GPD, by hand. Guided by Eq. (19), we would suggest for e, e.g., the ansatz: 

^ _ qE{y,t) N{b,p,a) ^(l-y)^ + y^-y(l-y)j [(i_^)2_^2]fe 
e(2/, z, t) - — — -— — - - ; ^ - - — — — p 7X9};^] • (35) 



^1 - [(1 -y)^ + y^-yii-y)- ((i - y? - " " ^)''^' 

The Regge improved PDF analog, generically parameterized by its large and small x behavior 



where qE{x,t) is normalized at t = to the anomalous magnetic moment k. The normalization 

N(h Pa)= r(fe + |)r(2-^(o) + /3)r(2-2P-^(t) + /3) 

^ ' V?r(6+l)r(2-2P-a(0)+/5)r(2-a(t)+/3) ^ ^ 

is introduced in such a way that the t-dependence in the form factor roughly factorizes as 1/(1 — 
a{t)), resulting from the Regge behavior, times an impact form factor, behaving as {—t)~^ for 
t —00 . The parameter b adjusts the power behavior of E{C,,^,t) ~ (1 — ^Y^^ at large ^ and 
fixed t. The parameters should be fixed from fitting of observables; their generic values read in 



11 



accordance with Regge phenomenology, counting rules [2, 3, 32, 53], and the spectator model: 

~ 0.5 + 0.9tGeV-^ /3 ~ 5, P ~ 2, 6-1, A~0.8GeV, m~0.4GeV. 

The model features, we spelled out here in momentum fraction representation, are also imple- 
mented in the Mellin space GPD ansatz [14, 15]. We emphasize that the ansatz (35) is still not 
optimal, since a flexible adjustment of the resulting magnitude for the amplitude at given t, in 
particular at small ^, is not incorporated and positivity constraints are no more automatically 
satisfied. A more detailed discussion of improved GPD ansatze should be given somewhere else. 
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FKZ 06 BO 103, by the International Cooperation Program of the KICOS (Korea Foundation for 
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